Manipulating the Tomonaga-Luttinger exponent by electric field modulation 
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We establish a theoretical framework for artificial control of the power-law singularities in 
Tomonaga-Luttinger liquid states. The exponent governing the power-law behaviors is found to 
increase significantly with an increase in the amplitude of the periodic electric field modulation 
applied externally to the system. This field-induced shift in the exponent indicates the tunability 
of the transport properties of quasi-one-dimensional electron systems. 
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I. INTRODUCTION 

Interacting electrons in one-dimensional (ID) met- 
als constitute a highly collective state of matter: the 
Tomonaga-Luttinger liquid (TLL) stated— The collec- 
tive nature of the TLL states is what distinguishes them 
from their higher-dimensional counterparts. Interact- 
ing electrons in two or three dimensions form a Fermi 
liquid, wherein the only effects of interaction are the 
modification of their effective mass and the possibility 
of being scattered. In one dimension, however, even 
the slightest correlation between electron motions has a 
dramatic effect, leading to distinctive features that can- 
not be explained by the Fermi liquid theory. To date, 
physical consequences of the TLL states have been ex- 
perimentally observed in various systems, including car- 
bon nanotubesj^"— semiconducting quantum wiresj^r— 
quasi- ID organic conductors quantum Hall edge 
states^ and other materials having highly anisotropic 
conductivity J^r— From the theoretical viewpoint, a more 
general TLL theory with a nonlinear dispersion^ as well 
as a novel wave-packed dynamics through Y-shaped TLL 
junctions^ have been recently suggested. 

A hallmark of TLL states is a pseudogap in the one- 
particle density of states v{e) at the Fermi energy £p. 
Injection of an additional electron in the TLL ground 
state disrupts the pre-existing correlation, thus requir- 
ing excitation of an infinite number of collective modes. 
This results in a power-law singularity of the form v{e) cx 
\e — £_f|", where a(> 0) is called the TLL exponent. 
The same power-law arises in the case of a differential 
tunneling current^ dl/dV cx \V\ a at high bias voltages 
(eV S> fcfiT) and a temperature-dependent conductance 
G(T) cx T a at low voltages (eV < k B T), although a 
may change due to environment effects.— These power- 
law behaviors are in strong contrast with the behavior of 
Fermi liquids; in the case of Fermi liquids, v{e) close to 



£p and dl/dV become constant. 



The TLL exponent a is nonuniversal; it is dependent 
on the interaction strength^ the geometric shape of the 
system^ and the position of tunneling^ In fact, dif- 
ferent values of a were obtained in carbon nanotubc ex- 
periments when each electron tunneled into the end or 
bulk of the system^£i2I A further non-trivial shift in a 
was suggested in multiwalled nanotubes, where a varies 
in a continuous manner under the application of a high 
transverse magnetic fields Continuous variation in a 
was also found in a nuclear magnetic resonance study 
of CuBr 4 (C5Hi 2 N)2 crystals;— in this case, an exter- 
nal magnetic field acted as the chemical potential. Such 
field-induced variations in a can be exploited for achiev- 
ing artificial control of transport properties in quasi-lD 
conductors, which would play a fundamental role in the 
development of next-generation quantum devices. 



In this paper, we propose a theoretical framework for 
TLL exponent manipulation based on electric field mod- 
ulation. An analytical expression of the exponent a for 
quasi-lD conductors subjected to a stepwise periodic po- 
tential is established in terms of the potential amplitude 
and period. Under feasible physical conditions, a in- 
creases significantly with an increase in the potential am- 
plitude; this indicates that it is possible to tune the quan- 
tum transport properties of quasi-lD systems by manip- 
ulating a. For a concise description, we focus our at- 
tention on ID spinless fcrmion systems, considering that 
the effects of spin degree of freedom requires no substan- 
tial revision of the present conclusion. This issue will be 
revisited in Sec. IV. 
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II. ONE-PARTICLE STATES IN ID PERIODIC 

SYSTEMS 

A. Eigenenergy analysis 

We consider a quasi- ID electron system having a thin 
cylindrical shape with length L, which is subjected to a 
periodic external potential field. The cylinder radius d 
is so small (d <C L) that all electrons reside in the low- 
est subband, Xj of the transverse motion. Single-particle 
wavefunctions thus have the form 



*(z,r ± ) = x{ri_)tp{z), 



(1) 



where the two-dimensional vector r±_ = (x, y) spans the 
circular cross section. The axial component ip(z) obeys 
the Schrbdingcr equation with effective mass m* , 



h 2 d 2 ^ 



u(z)ip(z) = eip(z), 



2m* dz 2 

where u is a stepwise periodic potential given by 

u(z) - [ °> ( ° " z < a) 

[ > ~ \ u > 0, (-6 < z < 0) 



(2) 



(3) 



and u(z) — u(z + a + b). The periodicity of u(z) implies 
that the eigenstates of Eq. ([2]) are represented by the 
Bloch function 



1 



(j>k{z)e 



Ikz 



b k {z + a + b)^(j) k {z), (4) 



where \tj}(z)\ 2 dz = 1 and J^ +b \<p(z)\ 2 dz = 1 with 
L = N(a + b). Hence, all eigenstates are labeled by the 
index k. It is noteworthy that the periodic potentials sim- 
ilar to the above can be realized by introducing geometric 
curvature (instead of field modulation) to quati-lD sys- 
tems, and such the curved systems may show non-trivial 
quantum transport ! 34 ' 35 

The dispersion relation of the system given by Eq. @ 
is (see Appendix A) 



/(e) = cos [k(a + b)} 
Here, the function /(e) is given by 



(5) 
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FIG. 1: (a) Variation of the curve /(e) with an increase in the 
potential amplitude from uo = 0.2 (solid) to 0.5 (dashed), 1.0 
(dashed-dotted) , and 2.0 (dotted) in units of Ry = e 2 /(ean) 
wherein as = h 2 e/(m*e 2 ). The parameters a = 5.0 and 
b = 1.0 in units of »b are fixed, (b) Slow decay of ef with 
increasing uq; see text for the definition of ef- 



radius ob = h 2 e/(m*e 2 ) for length and units of the ef- 
fective Rydberg Ry = e 2 /(eas) for energy, where e is the 
dielectric constant of the wire. For GaAs-bascd quantum 
wires, for instance, we have <2b ~ 10 nm and Ry ~ 50 
meV^ 

Figure QJa) shows the plot of /(e) corresponding to 
the lowest energy band of the allowed e. The poten- 
tial amplitude uq ranges from 0.2 to 2.0, and fixed pa- 
rameters a = 5.0 and b = 1.0 are used. With increas- 
ing Mo, the band width shrinks monotonically and the 
ground-state energy e gl -d (i.e., the specific e that gives 
/(e) — 1 in the plot) shifts to the right. Thereafter, 
we set the Fermi energy ep (measured from e gr d) such 
that the electron density n satisfies nae = 0.3, which 
gave r s = (2710b) 1 = 1-67. This value is much lower 
than the critical value of r s , viz., r s = 36, at which 
the Wigner crystal transition takes placed; it is also 
lower than r s ~ 2.2, above which the 4k f correlation 
in the charge density distribution was suggested^ Since 
n = (2/7r) Jq F dk, we can evaluate ep for a given uq from 
the relation /(ef) = cos[^n(a + b)]; see Eq. ([5]). Figure 
OJb) presents the M -dependencies of ep. It decays with 
increasing uo but remains of the order of 10 -2 Ry for all 
Mo under consideration. 



no 



rjo 

111 



sin(CaJ?o)sinh(Ch?7i) 



+ cos (CaVo) cosh (CbVi) > 



(G) 



where 770 = \J~£, r\\ = ^/uq — e, Ca = 2m*a 2 /H 2 , and 
C 2 = 2m*b 2 /h 2 . According to Eq. /(e) must fall in 
the range —1 to 1 for e to be physically relevant; i.e., only 
e's that satisfy the condition |/(e)| < 1 are allowed to 
be the eigenenergies of the system. Note that Eq. ((5|) re- 
duces to the trivial relation e = h 2 k 2 / (2m*) when w = 0, 
since Eq. ^ becomes f(e,u n = 0) = cos[(£ a + C,b)\fe\. 
Throughout this paper, we use units of the effective Bohr 



B. Field-induced shift in Fermi velocity 

The group velocity v g (e) = hr^de/dk of the system 
is evaluated by differentiating both sides of Eq. ([5]) with 
respect to k. The result is 



. + b y/l-f(ey 
h df/de 



(7) 



The explicit form of df/de is given in Eq. (|A5|) in Ap- 
pendix A. v g vanishes when e satisfies /(e) = ±1, i.e., at 
the lower and upper band edges; between the two edges, 
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FIG. 2: (a) Profiles of the group velocity v g (e) in units of 
vo = Ryae/ft. Each curve corresponds to the value of uo 
in the same manner as Fig. [T] A concave curve connecting 
the two end points is the trajectory of the Fermi velocity 
vf = v g (e = ef) obtained through the uo variation, (b) The 
Fermi velocity vf as a function of tto- 



where 

.94 = LV{k F ,k F ; q z =0), (10) 
g 2 = LV(k F ,-k F ; q z = 0), (11) 

9i = [V{k F , -k F ; q z = -2k F ) 

+V(-k F ,k F ; q z = 2k F )} (12) 

and 

V(h,k 2 ,q z ) = (k 1 +q z - k 2 -q z \V\k 2 ] h) . (13) 

The right side of (|T3|) is the matrix element of 
the Coulomb interaction between two-electron states 
(rj\kj) = ^{rj) given by Eq. ([T]), where q z describes the 
momentum transfer in the axial (z) direction. It should 
be noted that in formula (J9)), we can prove that g± = g 2 as 
a consequence of the symmetric property <pk{z) = (t>*_ k {z) 
of the function (j>k(z) introduced in Eq. See Eq. (f2"2")l 
for the proof. 



v„ takes the maximum value. The only exception oc- 



curs in the limit of uq = 0, in which v g (e) = ^J2e/m* 
increases monotonically with e. 

Figure EJa) shows the e-dependence of v g in units of 
vo = Ryae/fi. uq varies from 0.2 to 2.0 in the same 
manner as in Fig. [1] With increasing uq, the maximum 
value of v g decreases and the peak position shifts to a 
lower e. A concave curve connecting the two end points 
is the trajectory of the Fermi velocity v F = v g (e = e F ) 
obtained through the uo variation. The value of v F de- 
creases monotonically as uo increases, and it then con- 
verges to the origin in the limit of uq — > oo. Figure [5Jb) 
shows the plot of v F vs. uq. Note that v F /vo ~ 0.1 
corresponds to v F ~ 10 5 m/s if we employ the material 
constants of GaAs.— 



III. BOSONIZATION OF ID PERIODIC 

SYSTEMS 

A. TLL exponent in ID periodic system 

Here we discuss the effect of field modulation on the 
TLL states. An important indicator of TLL state realiza- 
tion is a power-law singularity of the one-particle density 
of states v{e) near e F represented by£ 



v{e) oc |e — e F \ a , a 



Here, a is the TLL exponent, and the parameter K, is 
defined as (see Appendix B) 



K 



2-nhvp + (g 4 - g 2 + gi) 
2ttHv f + (gi + g 2 - gi) 



1/2 



(9) 



B. Fourier representation of Coulomb interaction 

Computation of the TLL exponent a requires evalu- 
ations of V(ki,k 2 ,q z ) at the specific values of k\,k 2 ,q z 
given above. The explicit form of V(k\ ,k 2 ,q z ) is obtained 
by assuming the screened Coulomb potential 



ri,rj\V\rj,ri) = V{r. l:j ) = 



e -~i'. 3 i 



(14) 



where c c = e 2 /(4-7re) and k is the inverse of the screen- 
ing length. Equation (j!4[) has an alternative expression, 
given as 



(15) 



v/ (2tt)3 J ^|q| 2 + K 2 
From Eqs. (fTB")) and (fT5j) , we obtain (see also Eq. (jB6j) 

frtu i, ^ x [ dc t'z f rf< ?' 
V(k u k 2 ,q z ) = - \ — 



L J 2tt J (2tt) 2 
47rc e |F(ql)| 2 



G{q z , q' z , k 1 )G u (q z ,q' z , fc 2 ),(16) 



where A = a + b, q± = (q Xl q y ) is the transverse compo- 
nent of the three-dimensional wavevector q, and 

F(q' ± )=Jdr ± \ X (r ± )\ 2 e-^- r \ (17) 

G{q z ,q' z M) = Jdz^+gMl^Me-^-^m 

Gu(q z ,q' z ,k 2 ) = [ dz 2 ^_ q +z 2 )^ k2 (z 2 y^-^ z \ (19) 



The subscript U in Gjj and in Eq. (JTHJ) indicates in- 
tegration within the unit cell domain z 2 G [0,A]. It is 
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FIG. 3: (a) The TLL exponent a as a function of the field 
modulation amplitude wo. The cylinder radius d of the quasi- 
1D system is taken to be d = 0.1 (dashed), d = 0.3 (dashed- 
dotted) and d = 1.0 (solid) in units of as- The pamameters 
a — 5.0, b = 1.0, n — 10~ 4 are fixed for all curves, (b) in- 
dependences of <?4(= (72) (thick curves) and gi (thin) under 
the numerical conditions same as in (a). 



natural to assume that x( r J-) nas a Gaussian form such 



1/2 



exp 



\r±[ 
<P 



(20) 



which gives F(q' ± ) = exp (— \q' ± \ 2 d 2 /2) . As a result, the 
integral term with respect to q' ± in Eq. (|16p is rewritten 
as 



-kl! 2 d 2 /2 



(q' z ) 2 + W ± \ 2 + K 2 Jo t+(q' 2 + ^)d 2 



-« e {e^y Ei [- ( q ' z 2 + K 2 ) 



(21) 



where Ei(x) is the exponential-integral function defined 
by Ei(x) = — f^' x t~ 1 e~ t dt. Finally, we obtain the ex- 
plicit form 



V(ki,k 2 ,q z ) 



A 



x e 



4nL 



dq' z G(q z ,q' z ,k 1 )G u (q z ,q' z ,k 2 ) 

(22) 



' 2 t 2\ .2 

<7 Z + k id 



by which we can compute 174, 32, 5i in Eqs. ([Il3| - (|T2)) for 
a given hp. 



fixed the pamameters a = 5.0, = 1.0, and na-Q = 10~ 4 
so as to satisfy the condition n <C kp in accord to the 
bosonization procedure^ 

The most important observation in Fig. [3Ja) is the 
monotonic increase in a with increasing uq. In partic- 
ular, the data of a for d = 1.0 exhibits a pronounced 
enhancement. Such the uo-driven shift in a is caused 
by a simple mechanism based on the behaviors of con- 
stituents vp, gi(— 52) and g\ contained in formula Q 
of K,. First, it follows from Fig. [3jb) that (74 and g\ 
are almost independent of the change in uq and satisfy 
the inequality 34 > g%. Hence, a decrease in Vf causes 
a decrease in K, defined by Eq. ©. Second, the above 
inequality implies K, < 1, and thus the decrease in K. 
results in an increase in a as understood from Eq. ([5]l. 
As a consequence, a can be raised by enhancing the field 
modulation amplitude uq. For u$ 10.0, a converges to 
a limiting value determined bya = (/C + /C _1 )/2 — 1 and 
^ = [.9i/(2.94 — 5i)] 1 1 since «f vanishes asymptotically. 
The upper limit of a is dependent of the geometric pa- 
rameters d, a, b] therefore, it is crucial to set appropriate 
values of the three parameters in order to obtain optimal 
controllability of the exponent a. 

It should be remarked that we have discussed spinlcss 
fcrmion systems; they are realized in spin-polarized ultra- 
cold fcrmionic gases^i and in spin-1/2 antiferromagnetic 
ladders ^2r— whereas many quasi- ID systems involve the 
effects of the spin degrees of freedom. When we take into 
account the spin degrees of freedom, Eqs. ^ and © can 
be rewritten because of the presence of SU(2) symmetry 



1 



1 

r 



2-Khvp + (ff 4 || + g 4± - ggj - ff2_L + giy) 
2-nhvF + (.g 4 || + 34± + g 2 \\ + ,92_l - gi\\) 



(23) 

1/2 
(24) 



Here, g^\ and gt± (i = 1, 2, 4) express the matrix elements 
of interaction between electrons with parallel spin and 
those with antiparallel spin, respectively. In the usual 
Coulomb interaction, the matrix elements do not depend 
on the spin degree of freedom; consequently, g^ = g^± 
is expected. Therefore, we anticipate that the sizeable 
shift in a demonstrated in Fig. [3] should be qualitatively 
correct even if we take into account the spin degree of 
freedom. Details including the actual calculation will be 
shown elsewhere. 



IV. RESULTS AND DISCUSSIONS 



SUMMARY 



We demonstrate below that the value of a can be tun- 
able artificially by imposing an appropriate magnitude of 
the periodic external field. Figure 03a) shows the field- 
induced change in the TLL exponent a for different val- 
ues of the cylinder radius: d = 0.1 (dashed), d = 0.3 
(dashed-dotted) and d = 1.0 (solid) in units of ob- We 



We have theoretically shown that the TLL exponent a 
can be artificially controlled by applying a perodic exter- 
nal electric field to the system. The bosonization proce- 
dure has been used to obtain the analytic expression of 
the exponent a as a function of the amplitude uq and pe- 
riod A of the external field modulation u(z). The result 
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indicates that a increases significantly with the poten- 
tial amplitude Uq, whose magnitude is within the realm 
of the existing experiments. The significant variation in 
a is attributed to the field-induced shift in the Fermi 
velocity vp of the single-particle state. The present re- 
sults indicate that the quantum transport properties of 
quasi- fD systems can be tuned by manipulating a. Ex- 
perimental confirmation will provide a novel approach to 
manipulating fD quantum systems. 
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Appendix A: Dispersion relation of ID periodic 

systems 

Here, we derive the dispersion relation of Eq. (O for 
a ID system subjected to the stepwise periodic potential 
given by Eq. ©. It follows from Eqs. © and © that 
ip{z) behaves as 

J Ae iKz + Be~ iKz , (0 < z < a) 

CeQ' + De-*, (-b<z<0) (A1) 

and ip(z + a + b) = ^(z)e ik{a+b \ where e = K 2 K 2 /{2m*) 
and uq — e = h 2 Q 2 / (2m*). The continuity conditions of 
ip and dip / dz at both z = and z = a lead to the matrix 
equation M v = 0, in which v = [A, B, C, D] T and 



M 



1 



-1 



f f 

iK -iK -Q Q 

a 1 1 a —u/p —pui 

iKa —iK/o —Qui/p Qpui 

oQb 



(A2) 



e""*, p = e^", and lu = e ik i a +t>). Solving the 
secular equation detiVf = 0, we can obtain the desired 
result of the dispersion relation ([5]). 

In general, the determinant detA<f of a square ma- 
trix M consisting of n 2 elements {a^} is represented 
by detM = Y^=i(~ lT +j a ijC>ij with an arbitrary fixed 
value of j, where CV, is the determinant of the (n — 1) x 
(n— f ) matrix that is obtained from M by erasing its zth 
row and jth column (see Ref4£ for the proof). Hence, 



for Q ^ 0, the determinant of matrix M given by Eq. 
[1 is 

detM 



QqoJ 

iK 
a 



I 


1 


-1 




1 


I 


-1 


iK 


-iK 


-Q 




iK 


-iK 


Q 


a 


1/a 


-u/p 


P 


a 


1/a 


-pu> 


I 


-I 


-1 




f 


-1 


-1 


iK 


-Q 


Q 


— iKa 


-iK 


-Q 


Q 


a 


-u/p 


—puj 




1/a 


-u/p 


—puj 



MKQ (f + uj 2 - 2uj cos Ka cosh Qb) 
+ii (K 2 - Q 2 ) uj sin Ka sinh Qb, 



and thus detiVf = reads as 



cos Ka cosh Qb - 



K 2 -Q 2 
2KQ 



(A3) 



sin Ka sinh Qb. 



(A4) 

This result is equivalent to the set of Eqs. <j5j> and ((6]). In 
addition, the derivative df /de is obtained by differenti- 
ating the right side of Eq. ([6]) with respect to e. Straight- 
forward calculation yields 



dj_ 

de 



% 



4?7i 



i 



hi 



Ka 

Cb ' 
20, _ 

Ca " 

vi 



sin(C a r?o)cosh(C b r?i) 
cos(C a 7/ )sinh(Ch??i) 
2 | sin(Ca??o)sinh(Cf ) 77i).(A5) 



Specifically, when Q = (i.e., e = uq), the continuity 
conditions of tp and dip/dz at z = a yield Ka = ^7r and 
k(a + b) = £tt with an integer £. Given uq > 0, therefore, 
the eigenstate exactly at e = uq exists only if there is an 
integer £ that satisfies H 2 /(2m*) ■ (£ir/a) 2 = Uq. Other- 
wise, no eigenstate lies at e = u . 

Once an eigenenergy e and the corresponding k are 
obtained, we introduce them into the equation Mv = 
to evaluate the coefficients A, B, C, D of the eigenfunc- 
tion tl'(z) that belongs to e. The explicit forms of the 
coefficients are 



B 
A 

C 
A 



Q + iK e 



ik(a+b) —C 



iKa 



Q-iK e *Ha+b) e -Qb _ e - 

Q + iK i sin Ka 



iKo 



Q 



?ik(a-\-b) g — Qb ^—iKa ' 



(A6) 
(A7) 



and D/A = 1+ (B/A) - (C/A); A is uniquely determined 
by the normalization condition of ip(z). 



Appendix B: Proof of Eq. © 

The TLL parameter fC introduced in Eq. ((SJ) is deduced 
from the momentum representation of the two-electron 
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interaction Hamiltonian Hint for a quasi-lD periodic sys- 
tem. In the real-space representation, H; n t is given by 

Hint = \ f dm f dr 2 V(\m-r 2 \) 



x* t (r 1 )* t (r 2 )*(r 2 )*(r 1 ), (Bl) 



where 



*(n) = X] a fe' ! / ) fe(2i)x(r'i 



(B2) 



with a fermionic annihilation operator Cfc. Substituting 
Eq. (|B2|) into Eq. (fBl~|) . we have 



Hint = r X] ^fci,fe2;fc3,fe44 1 4 2 £ fe3Cfc 4I (B3) 



fci)fea,fc3,fc4 



and 



V kl ,k 2 :k 3 ,k 4 = J T2r± J dzi J dz 2 V (\r 1 -r 2 \) 

x r kl ( z iMM)^ 3 (z 2 )^ ki ( Zl ), (B4) 

in which the integration operator J T>rj_ = 
jdr\jdr\ \\ (r]_) | 2 |x (r^) | 2 acts on V r (|r 1 -r 2 |) 
since depends on r l ± . In the meantime, we derive 
an alternative form of Vfci,fc 2 ;fc3,fc4 (i.e., Eq. (|B6[) ) which 
allows to obtain the momentum representation of Hint 
(i.e., Eq. (jBllj) ) relevant to the formula of fC. 

Bearing in mind the periodicity of a system with period 
A = a + b, we replace the variables Z\,z 2 (0 < Zi < L, 
i = 1,2) in Eq. (JED with zi + /'A and z 2 + (Z' - Z)A 
(0 < 2^ < A, i = 1,2), respectively; the integers I and Z' 
are independent of each other and vary from to N — 1 . 
Applying the periodic boundary condition of ipk(z) — 
ipk(z + L) for an arbitrary k, we obtain 



Vi 



k 1 ,k 2 ;k 3 ,k 4l 



1 

iV 5 



EE 

l V 



Vr, 



dz\ I dz 2 



x V (\m-r 2 +lXe z \)^(zi,z 2 ) 

X e ~ i ( fe i+ fc 2-fc3-fc4)C ^ e i(k2-k 3 )lX (B5) 

where (j){zi,z 2 ) = ^(^l)^ (22)^3 ( 2 2)^fe*(«i) and /„ 
symbolizes the integration within the unit cell z = [0, A]; 
e z is a unit vector in the axial direction. Summation 
over I' in Eq. (|B5|) yields a term NSk 1 +k 2 ,k 3 +k4 ( we con- 
sider only normal scattering contributions, omitting the 
Umklapp scattering processes). We further assume a 



momentum conservation that requires relabeling of the 
set (fcx, k 2 , &3, fej) by (ki + q,k 2 — q,k 2 ,k\). Introduc- 
ing the assumption into Eq. (|B5j) and performing sum- 
mation with respect to Z, we find that V(ki,k 2 ; q) = 
Vki+qM-q-MM is given by 

V(k!,k 2] q) = ^Jvr ± Jd Zl J dz 2 V (\m-r 2 \) 
x^ i+ ,(z 1 )fe(z 1 ) e -^4_,(z 2 )fe(z 2 )e'«(B6) 

In Eq. (|B6|) , the subscript U originally attached to J v dz\ 
vanishes so that z\ ranges from to L. 

We are ready to evaluate the TLL parameter JC. From 
Eq. dB6j) and Eq. dB3j), it follows that 



^int = \ ^2V(ki,k 2 ; q)c{ 1+q c\ 2 _ q 



(B7) 



ki,k2,q 



Substituting ki = ptkp + n into Eq. (|B7[) under the as- 
sumption that of I At I <C kp and p = ±1, we get 

Hint = ^E E [v(pk F ,pk F ; 0)4,«+g4,«'-g^.«' a P." 

p k,k' ,q 

+ V(pk F ,-pk F ; -2pkF)cl lK+q cL PtK ,_ q c PtK >C- p , K 
+ V{pk F ,~pk F ; 0)^ iK+q c!_ PtK ,_ q C- PjK >dp lK , (B8) 

where c p . K = c p k F + K - Since the time reversal symmetry 
of the system assures that Vk^-MM = V-k lt -k 2 --k 3 -k^ 
we have V(k\,k 2 \ q) = V(— k\ — q,—k 2 + q; q) that 
implies 



V(k F ,k F ; 0) = V(-k F ,-k F ; 0), 
V(k F ,-k F ; 0) = V(—k F , k F \ 0). 

As a consequence, we obtain the result 



%l iL \ " \ " at at - 

rt int — 2^ Z-f 2-^1 C P,K+q C p,n' -q C P,K' C p,K 
V K 7 n',q 

92 - .9i 



(B9) 
(B10) 



E E 4,K+g a -p,K'-?C_p, K 'Cp, K ,(Bll) 
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where 34, _g 2 , and #1 are defined by Eqs.(|T0 |) -(|T2 |) . 
Through the diagonalization of Hint in Eq. (|B11[) . we fi- 
nally obtain the TLL parameter JC given in Eq. ([9]) ; this 
parameter depends on the two coefficients g 4 and g 2 — g\ 
prior to the summations in Eq. (|Blip . 
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